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Abstract 

A comprehensive theoretical model for the bulk manganite system Lai_ x (Ca,Sr) x Mn03 is pre- 
sented. The model includes local and cooperative Jahn- Teller distortions and the on-site Coulomb 
and exchange interaction. The model is is solved in the single-site dynamical mean field approx- 
imation using a solver based on the semiclassical approximation. The model semi-quantitatively 
reproduces the observed phase diagram for the doping < x < 0.5 and implies that the manganites 
are in the strong coupling region but close to Mott insulator /metal phase boundary. The results 
establish a formalism for use in a broader range of calculations, for example on heterostructures. 

PACS numbers: 71.10-w,71.30.+h,75.10.-b 
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I. INTRODUCTION 



LaMnOs crystallizes in a structure closely related the basic ABO3 perovskite form. As 
the temperature is varied it undergoes orbital ordering and antiferromagnetic transitions. 
Replacing some of the La by divalent alkali ions such as Ca yields an even wider range of 
phenomena, including charge ordering, ferromagnetism, and colossal magnetoresistance [? ? 
? ]. The phase diagram is summarized in Fig([T]). While the manganites have been studied 
for many years, and much of the physics has been understood, there is as yet no consensus in 
the literature on a model which is rich enough to account for all the physics, includes all of 
important interactions, and can be solved to predict (or at least explicate) new phenomena 
such as those occurring in systems such as heterostructures. In this paper we develop such a 
model and explore its properties. Our results place the materials in the strong/intermediate 
coupling regime. 
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FIG. 1: The experimental phase diagram as a function of doping x and temperature T. PO: 
paraorbital; PM: paramagnetic; FM: ferromagnetic; OO: orbitally ordered; AAF: A-type antifer- 
romagnetic. See text for the descriptions of the phases. 

The phase diagram shown in Fig([T]) includes two magnetic orders, A-type antiferromag- 
netic (A-AFM) and ferromagnetic (FM) states, corresponding to the Mn spin arrangements 
shown in FigQ. The A-AFM structure consists of ferromagnetic planes antiferromagneti- 
cally coupled. With our coordinate choice, each FM plane is spanned by x and y while the 
remaining direction is z. The orbital order (OO) in this context refers to a particular distor- 
tion arrangement where the oxygen octahedra have in-plane staggered (x — y plane) Q x JT 
distortions plus a uniform — Q z distortion (the minus sign represents the octahedron shrink- 
ing in z while expanding in x — y directions) (FigQ). The metal/insulator phase boundary 
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FIG. 2: (Color online) Illustrations of FM, A-AFM, and orbital orders. The filled and open circles 
represent manganese and oxygens respectively. The light arrows stand for the core-spin orientation 
at each plane. Upper left panel: illustration of two magnetic orders. For FM, the spins at different 
planes are aligned in the same direction while for A-AFM, the spins at adjacent planes arranged 
oppositely. Upper right panel: illustration of the in-plane staggered Q x order is shown in the up- 
right corner. The double-arrows represent the long O-Mn-0 distance caused by the Q x distortion. 
Lower panel: illustration of the uniform —Q z order: the system uniformly shrinks in the z direction 
while it expands in x — y. 

is determined from the DC resistivity. The definition of metal/insulator is ambiguous. Here 
we adopt the definition that the system is metallic/insulating at a given temperature if the 
temperature derivative of resistivity is positive/negative. 

In this paper, we present a model which captures all of the physics discussed above and 
solve it by the single-site dynamical mean field theory (DMFT)[? ]. There are two main 
purposes for this study. First, although basic understanding for exhibited phases at a given 
doping is known, it is important to determine the extent to which the general model with a 
fixed set of parameters matches the observed phase diagram. Second, we wish to apply 
this theory to understand the behavior of the recently synthesized manganite superlattices 
[? ? ]. 

Solving a theoretical model ordinarily requires approximations. Here we use the single 
site DMFT[? ]. This approximation requires as an intermediate step the solution of a 
quantum impurity model. In this paper we solve the impurity model using a generalization 
of the semiclassical approximation (SCA) [? ]. We generalize it to the 2-band case and 
develop a formalism for incorporating the cooperative Jahn- Teller ( JT) effect into the single- 
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site DMFT. We semiquantitatively reproduce the observed phase diagram for < x < 
0.5 and identify the sources of the observed phases. Our calculation yields three main 
results. First, our calculation suggests the problem is in the "strong/intermediate" coupling 
regime in the sense that under the single-site DMFT approximation the local interaction 
strength is slightly larger than the critical value needed to drive a metal-insulator transition. 
Consequently the system is very sensitive to the mechanisms governing the bandwidth such 
as magnetic order and details of crystal structure. Second, the cooperative Jahn- Teller effect 
is the main source accounting for the observed high orbital ordering temperature. Finally, 
our calculation confirms that when the doping is increased to the colossal magneto-resistance 
(CMR) region x ~ 0.3 (N ~ 0.7), the double-exchange (DE) mechanism becomes dominant. 

The rest of the paper is organized as follows. We first present the model and the in- 
teractions included. After providing key steps for our approximation, we show how we fit 
parameters and then present the results. Discussion concerning the inadequacies of our 
model/approximation and differences between the calculation and the experiments is given 
in Section VII. Section VIII is a conclusion. In the appendices we examine the validity of 
the SCA approximation and give details of the procedure we use to take the cooperative 
Jahn- Teller effect into account. 

II. MODEL HAMILTON! AN 

In this section we describe the interactions included in our model and define terms and 
notations which shall be used for the rest of the paper. 

Tight-binding: The band structure is described by a tight-binding model where only near- 
est neighbor hopping between e g orbitals is included. A justification for this approximation 
is given in [? ]. Two e g orbitals are labeled as |1) = \3z 2 — r 2 ), |2) = \x 2 — y 2 ). This implies 
a band Hamiltonian which may be written as 

Hband = ^2 € k,ab,a C l :a>a %b,a i 1 ) 
k,ab,cr 

e kaba = ~^( £ o^ + £ z^z + £xT x )ab where f 's are Pauli matrices, e is the unit matrix and 
eq = cosk x + cos k y + cosk z , e z = cosk z — ^(cosk x + cosk y ), and e x = ^{cosk x — cosk y ). 
a,b label orbitals, i,j sites, and a spins. We emphasize that what denoted here as two e g 
orbitals are actually the anti-bonding combination of Mn 3d and its neighboring oxygen 
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2p(cr bond) states [? ]. 

On-site e-e: For the on-site interaction within e g orbitals, we use the Goodenough- 
Kanamori-Slater approximation in which the form of interaction is the same as in the free 
atom. Two independent parameters conventionally denoted as U and J are required to 
specify this interaction. It is generally accepted [? ] that the charging energy U may be 
strongly renormalized by solid state effects whereas the inter-orbital exchange energy J is 
less affected. The e-e interaction within the e g multiplet is 

H EE = /.{U ~ J) n l,v n 2,v' + U 22 n i,t n hl + J{ C 1,T C 1,J. C 2,1 C 2,T + h ' C ) ~ 2 *^1 ' ^2 (2) 
a,a> i=l,2 

where £1(2) = ff af 3c\^ a ci(2),p- The J( cj^c^j^j^-r + h.c.) term is referred to as the pair 
hopping and the — 275*1 ■ 5*2 term is the exchange. 

Hund's coupling: The coupling between Mn e g and Mn t 2s electrons is approximated by 
H Hund in which three t 2g electrons are treated as an electrically inert "core spin" of magnitude 
S = (3/2)h. We shall further approximate the core spin as classical and normalize J# by 
taking |5| = 1, leading to 

HHund = — Jh Si ■ C i a (J a pCi t p (3) 
i 

where Jh > and |5| = 1. The minus sign ensures that the high spin state is energy-favored 
in accordance with the Hund's rule. 

Lattice elastic energy: For the lattice degree of freedom, we consider Mn motions in 
arbitrary directions^) and oxygen only along Mn-0 (a) bond direction (ui iX ) [? ? ] where 
5{ and illustrated in FigQ. The general lattice elastic energy in the harmonic 

approximation is 



Hiat = EK^ a - <? + W - <-J 2 ] (4) 
2K M _ ^ 

1, a 

+ \ E Eab m^ + \ E Dab m u -k ( 5 ) 



where 1/Km-o is spring constant between neighboring manganese, oxygen, while E ab (k) and 
D ab (k) are general Mn-Mn, 0-0 couplings in k space. In the specific numerical calculations 
presented here, we set E ab (k) = and D ab (k) = fc 4 ^ 5 ab sin 2 (A; a /2), but effects arising from 
a more general interaction are discussed. With our convention, u, 5 and lattice constants K 
all have dimension of energy. 
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FIG. 3: (a) The ideal cubic perovskite structure for LaMnOs and the lattice degrees of freedom 
considered here: Mn can move in arbitrary direction 5i while oxygen ions only move along the 
Mn-0 bond direction u^ x r y ^ z y (b) Sketch of the three octahedral distortion modes considered here: 
Qo breathing mode, Q x and Q z Jahn- Teller modes. 



Electron- Lattice coupling: The breathing (Qo) and Jahn- Teller (JT) (Q x , Q z ) modes at 
site i are denned by: 

Qifl = —^(vi, x + v it y + v i>z ) 



Qi,x x ^i>y^ 



Qi, z = —f=(-v itX - v i)y + 2v i)Z ) (6) 
where v^ a = u^ a — Ui- a , a . The e g orbitals couple to these three modes as 

HjT = -J2(^ T ab + QiXb)4,a^b (7) 
i,a,b 

H B = -PQi,o(ni - (n)) (8) 

where (3 is dimensionless and positive. We will take (3 = 1 which simplifies the discussion of 
cooperative Jahn- Teller effect. In this paper we treat Hb by a mean field approximation, so 
it is only important when the charge distribution is not uniform, like in the heterostructures 
or in the charge-ordered phase. In our definition, positive Q z stands for the distortion where 
the octahedron expands in the z direction while shrinking uniformly in x — y with fixed 
volume. The minus sign in Eq([7]) means that the positive Q z favors the occupancy of the 
\3z 2 — r 2 ) state. This sign choice is justified because positive Q z increases the lattice constant 
in z direction and consequently reduces \t p d\ and E anti _ bonding , increasing the occupation in 



the anti-bonding band which is mainly composed of Mn \3z 2 — r 2 ). A similar consideration 
leads to the minus sign in Eq(jHJ) (positive Qo means a volume expansion of the octahedron). 

Cubic term in lattice energy: An anharmonic cubic term [? ] in lattice energy is also 
included. 

Hcmc = -A(3Ql - Ql x Qi, z ) (9) 

where A in our convention has the dimension [E] -2 . Note this is the only cubic combination 
satisfying the lattice cubic symmetry. With the minus sign, positive A is required to produce 
the observed distortions for LaMnOs. 

G-type AF coupling: There is an isotropic nearest neighbor AF coupling (G-type) between 
t 2g spins Si. 

Haf = Jaf Si • Si + h (10) 

i,n 

with positive Jaf- This coupling arises from the super-exchange mechanism (virtual hopping 
in t2 g channels) and experimentally shows in the G-type AF order exhibited in CaMnOs [? 
]. The main effect of this term is to reduce the magnetic transition temperature. 
The total Hamiltonian is then 

Htot = Hband + H EE + Huund + Hi at + Hjt + Hb + Hcubic + Haf (H) 

III. METHOD 

We use the single-site Dynamical Mean Field Theory (DMFT) with the semiclassical 
approximation (SCA) to solve this two-orbital problem [? ]. In the DMFT approximation 
one replaces the full lattice self energy E(u;,p) by a local (momentum-independent) quantity 
E(a>) which is determined from the solution of an auxiliary problem (quantum impurity 
model) plus a self consistency condition. The multiplicity of orbitals and interactions means 
that the impurity model is not easy to solve. We use a Hubbard-Strotonovich transformation 
proposed by Sakai [? ] and the semiclassical approximation. To evaluate the frequency sum 
we use a procedure recently introduced by Monien [? ]. 

We also mention two simplifications here. First we do not take into account the Coulomb 
potential produced by the random distribution of cations - the only effect of replacing 
some La by divalent elements is to reduce the e g electron population. However due to 
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the screening effect from conduction electrons, we believe this simplification is not crucial. 
Second, we restrict our calculation to charge- uniform states therefore we cannot obtain the 
charge ordered phase which may be energy favored around half doping. 

The rest of this section is organized as follows. Two key ingredients will be discussed: first 
we show how we encode the cooperative Jahn- Teller effect in the local impurity problem; 
second we give some detailed formalism about SCA in this 2-orbital problem, especially how 
we decompose the quartic interaction and what simplifications we have made. Then we 
discuss what measurements we used to fit parameters. 

A. Cooperative Jahn- Teller 

The local octahedral distortions (Qo,Q x ,Q z ) at different sites are not independent - 
distortion at one site inevitably causes distortion at the neighboring sites so that some 
global configurations of the lattice distortions are energy favored. This is the cooperative 
Jahn- Teller effect [? ] which correlates the octahedral distortions at different sites. Here we 
include this inter-site effect into the single-site DMFT by integrating out all of distortion 
fields except for those involving the variable v at the site of interest. The detailed calculation 
is given in the appendix A and the resulting local effective potential is 

V eff (Q ,Q x ,Q z ) = ^ + eF-Q (12) 

where K is an effective spring constant, eF represents the force exerted on the distortions at 
one site by static (mean field) distortions on the other site. Here F measures the amplitude 
of the long ranged order and e gives the strength of the cooperative Jahn- Teller coupling. 

B. On-site e-e 

The key step in our solution of the impurity model is to rewrite the quartic interaction 
into sums of complete squares so the continuous Hubbard-Strotonovich transformations can 
be applied. Using the decomposition proposed by Sakai [? ], we define f a = c\ a C2 a + cl a ci a , 
n = m +n 2 , q = n\ — n 2 , s= (n iiT — m,j.) + (n 2 ,f — n 2 ,i), and d = (m,f — ni,j.) - (n 2 ,t —^2,4.), 
and re-express Eq (J2J) as 

Hee = U n - J -{ f] - hf + ^ - ^ - ^ - ^ (13) 
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with U = 0,U n = (3U - 5 J)/4, U q = (U - 7J)/4, U s = (U + J)/4, and U d = (U - 3 J)/4. 
Due to the fermionic identity = nj i(7 (i = 1,2 <r =T, I), those coefficients are not unique. 
For example, C/ = J/2, C/ n = (3C/ - 6J)/4, U q = (U - 6J)/4, [/, = ([/ + 2J)/4, and 
Ud = (C/ — 2J)/4 is another legitimate set of choice. If the impurity problem is solved 
exactly, these two choices lead to the same result, but if approximate methods are used this 
needs not to be the case. However in the current study, the coefficients will be determined 
by fitting to data so this ambiguity is not important. 



C. The Impurity Problem 

The impurity problem is then described by the effective action S = So + Si nt where 

So = ~ J drdr'a%{r - t')4(t) C/3j (t') (14) 

and S int = J dTH EE {r). The partition function is Z imp = J dQ dS D[c* c] e~ 5 . Applying the 
Hubbard-Stratonovich (HS) transformations [? ] to decouple Hee, one arrives 

- 2TT s ^ T) + 2U d 
+ Jdr {-i<t> n {r)n{r) + 0/(r)(/ T (r) - f^r)) + [(fi q (r)q(r) + (fi s (r)s(r) + &(r)d(r)]) 

+ J dr(Q- f ab 5 a/3 + J H S ■ 8 ab a a ^j 4 a (r)c^(r) (15) 

To maintain the symmetries of the local interaction (SU(2) for spin, U(l) for orbital), we 
generalize the scalars (fi s , (fi q , to vectors (fi s (= (<j) StX , (fi s ,y, <fi s ,z), 3 components) and (fi q (= 
(<f>q,z, <Pq,x), 2 components) and average over their directions[? ]. After expressing the Si nt 
in frequency space, two simplifications are made. First, only zero frequency component for 
each HS field ((fi(iu>o) = (fi) is kept and second, saddle point approximations are applied to 
(fif, (fid, and Q fields, i.e. (fif = (fid = Qo = 0. Different methods have been proposed for 
handling the i(fi n field [? ? ]. For the 2-band model studied here, we found the method 
of Ref[? ] effectively enhances the local orbital moment as the doping increases which is 
opposite the observation while the method in Ref[? ] is free from this trouble, therefore we 
follow Ref[? ] and take i(fi n = 0. After integrating out the fermionic degree of freedom and 
combining the lattice effect, one gets 

Veff = + J- + J- ) + A(3(Q Z ) 2 - {Q X ) 2 )Q Z + eF ■ Q - T £ Tr log A(iu n ) (16) 



with 

A = a+(Q + (f q )-T+(J H S + s )-a (17) 

where A is a 4 x 4 matrix and a the Weiss function. The A(3(Q Z ) 2 — (Q X ) 2 )Q Z term comes 
from the simple mean field approximation of the cubic lattice energy. 

D. Parameters and Fitting 

The discussion above indicates that there are seven parameters to be determined: the 
hopping t, effective local JT coupling Uq, effective magnetic coupling U s , Hund's coupling 
Jh, core-spin AF coupling Jaf, cooperative JT coefficient e, and anharmonic lattice energy 
A. The first six have the dimension of energy (E) and will be measured in units of the 
hopping t while the last one has dimension 1/E 2 . 

The hopping strength t has been determined from the band structure calculation for 
the experimental observed structure of LaMn0 3 to be roughly 0.5eV (~5000 K) [? ]. By 
contrast, the value appropriate to the ideal perovskite structure is roughly 0.65eV. The 
difference is due mainly to the effect of the GdFeC>3 rotation. For the series La!_ x Ca x Mn03 
the distortion depends weakly on the doping x (less than 10%). For Lai_. E Sr. E Mn03 the 
rotation angle is more x-dependent (up to 30%). Therefore calculations in which t is taken 
to be independent of x may be appropriate for the Ca series but are unlikely to be adequate 
for the Sr series. We focus here on the Ca series. We note that the value t = 0.5eV is in 
good agreement with the spectral weight inferred from the optical conductivity experiments 
on the ferromagnetic phase of La .7Cao.3Mn03 [? ]. The Jaf is estimated from the Neel 
temperature of CaMnOs, roughly 110A ~ O.OleV [? ]. The T/v of the Heisenberg model 
H = Jaf\Si • Sj with \S\ = 3/2 obtained from simple mean field is ^Jaf, from which 
we estimate J af ~ 4.5meV~ 0.009t. The super-exchange argument (from virtual hopping 
of t 2g electrons) [? ] also leads to the same estimate. We further found that within the mean 
field approximation, J af term only acts to reduce to magnetic transition temperatures 
by 5-10% but does not result in any new magnetic order. Including this term however 
substantially increases the calculation time, so we typically set Jaf = to accelerate the 
converging processes. The anharmonic lattice energy is taken to be A = 0.006t~ 2 so that our 
calculation reproduces the observed JT distortion around room temperature for LaMn0 3 [? 
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The remaining parameters J#, U s , Uq, e are fitted by comparing the calculated and 
observed optical conductivity in LaMnCV Generally the <j(uj) contains peaks corresponding 
to local excitation energies of the system. If we treat the hopping t as a small perturbation, 
then peaks in the optical conductivity are roughly the energy differences between the excited 
states and the ground state of the local Hamiltonian. Since there are 4 states (2 spins x 
2 orbitals) at each site, we expect there are 3 main peaks in <j{oj) corresponding to three 
2-electron final states. The saddle point estimate from the local potential indicates these 
three peaks are located at 2Uq{\ + e) (correct spin, other orbital), 2(J# + U 3 ) (same orbital, 
antiparallel spin), and 2Uq(1 + e) + 2(Jh + U s ) (other orbital, antiparallel spin) which 
essentially agrees with our calculated results shown in Fig(jS]). The issue is discussed further 
in Ref[? ]. Experimentally there are two apparent peaks observed in LaMnOs [? j 
the lower one around 2eV (~ 4t) while the higher around 4eV (~ 8t). There are several 
minor structures around 5-6eV which we do not consider. Fitting the two main peaks in 
optical data suggests Uq{1 + e) ~ 4(t), J H + U s ~ 8(t). We choose J H = 2.8, U s = 1.4, 
Uq = 2.1, e = 0.05. We found that if Jh + U s is fixed the relative values of Jh and U s 
do not change the result much as long as Uq > U s . If Uq < U s , the orbital order is not 
stable against the magnetic order. Since at T = the SCA reduces to the simple mean field 
approximation where expectation values are determined by their saddle point values, the 
combination Uq(1 + €) uniquely determines the T = phase. However e has more significant 
effect on the non-zero temperature phase. We choose e so as to produce the observed T OD . 

Standard parametrization: Unless indicated otherwise the results we shall present later 
correspond to our standard parameters t = 0.5eV, Uq = 2. It, U s = lAt, Jh = 2.8t, Jaf = 0, 
e = 0.05, A = 0.006/t 2 . We remind the reader that this choice of t is only appropriate for Ca 
doped materials; in the Sr series the t changes with doping. All the temperature, frequency 
are measured in t; a simple conversion 0.5eV~ 5500K. 

IV. THE PHASE DIAGRAM 

The calculated phase diagram as a function of doping and temperature is shown in Fig(j4]) 
The results qualitatively reproduce the observed phase diagram (FigQ) in the sense that the 
relative positions of calculated magnetic/orbital phases are consistent with the experiments, 
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but the temperature scales are larger than observed. 
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FIG. 4: Calculated phase diagram as a function of doping x and Temperature. M and I stand for 
metallic and insulating phases. The dashed curve is the PM/FM phase boundary, computed using 
t = 0.6eV appropriate to Lao.6Sro.4Mn03. 

According to the doping, we divide the phase diagram into three regions - the undoped 
case (x = 0, LaMn0 3 ), the CMR regime (0.3 < x < 0.5), and the crossover regime (0 < x < 
0.3). In essence, the undoped case is cooperative JT dominated and the signature is the 
lattice distortions or equivalently the orbital order. On the other hand, the CMR regime is 
double-exchange (DE) dominated in which the system is FM/metallic at low temperature. 
In the crossover regime, both mechanisms play non-negligible roles to the system and we see 
that as doping increases, the cooperative JT effect decreases (T OQ decreases) while the DE 
mechanism gradually takes over (T c increases). In principle, we can extend our calculation 
to x > 0.5. However this region the effect of G-type AF coupling J AF starts to emerge (or 
both double-exchange and cooperative JT effect decreases) and a different self-consistency 
condition (G-AF) is required, so we shall leave it for future study. 

The remainder of the paper is organized as follows. We shall devote one section for the 
undoped case and one for the CMR and crossover region for more detailed discussions, then 
discuss the discrepancies between calculated and observed results. What important physics 
we are missing in our model/approximation and their effects to the current results will be 
stated. We also point out here that for the spectral functions we shall present, the fermi 
energy is at zero. Without further indications, p, a and A stand for resistivity, conductivity, 
and spectral function respectively. 
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V. THE UNDOPED CASE 

A. Overview 

Experimentally LaMnOs is insulating for all temperatures at least up to 800K [???], 
which is slightly greater than the orbital ordering temperature T ao . When the temperature is 
lowered, it first goes from PO/PM to OO/PM at T QO ~ 7801T(0.135*), then from OO/PM to 
OO/AFM at Taaf ~ 1401^(0.04*) where temperatures in Kelvin are experimentally derived 
while numbers in parentheses are from calculation. After a discussion about the nature of 
the insulating behavior, we examine the physical origins of the exhibited phases. 

B. The Transport and Excitation Spectrum 

in 

0.08 
_u 0.06 

'a 

1 004 
b 

0.02 



"0 2 4 6 8 10 12 14 11 "0 2 4 6 8 10 12 14 16 

Frequency (t) Frequency (t) 

(a) (b) 

FIG. 5: Optical conductivities for (a) T = 0.16*(> T 00 ) (dotted) T = 0.07 ~ 0.5T oo , OO/PM 
phase and (b) T = 0.02t ~ 0.5Taaf, OO/AAF phase. The heavy solid/dashed curves represent 
the in-plane/out-of-plane optical optical conductivities. To convert the frequency scale to physical 
units we note that the band theory indicates t = 0.5eV so * = 4 corresponds to 2eV. 

In this subsection we present our calculated optical conductivities for the stoichiometric 
end- member LaMnOs. We show that the calculated conductivities are in good agreement 
with experiment and that the agreement implies that LaMnOs is a Mott insulator; a more 
careful discussion regard to Mott insulator is given in Ref[? ]. 

To establish our claim we present in FigdS]) the optical conductivities for electric fields 
parallel to the x — y plane (solid line) and perpendicular to it (dashed line) at temperatures 
T = 0.16* (greater than T oa ), T = 0.07* (below T OQ , above the magnetic ordering tempera- 
ture Taaf), and T = 0.02* (roughly O.BTaaf)- The integrated optical conductivities upto 
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3eV qualitatively agree with experiments of LaMn03 [? ]. At T = 0.16t we see that the 
conductivity has two peaks at u ~ At and 8t, and a soft gap at u — 0. If we suppress 
the orbital order, forcing PO/PM solution down to lower temperature, we find that the low 
frequency conductivity decreases [? ]. When the temperature is lowered to T = 0.07t where 
the orbital order is well established; the peak positions remain essentially unchanged. An 
anisotropy produced by the orbital order appears and the gap at low frequency becomes 
sharper. 

As the temperature is further decreased into the A-type antiferromagnetic state the peaks 
sharpen and the anisotropy becomes more pronounced with an increase in a xx and a decrease 
in a zz for ui ~ At and the converse behavior in the u ~ 8t regime. This qualitative behavior 
was used by authors of Refs [? ] and [? ] to identify the lower feature as the transition 
to the maximal spin, orbitally disfavored final state and the higher feature as the transition 
to a lower spin, orbitally favored state. We make the same identification here and have 
adjusted the crucial parameters Uq and U s to place these peaks at the experimentally correct 
energies. Referring now to Fig(jSJ)(a) we see that for these parameters the correlations are 
already strong enough to produce an insulating state in the absence of the long ranged order 
which is one characteristic of Mott insulator, although the "soft" nature of the gap places 
the materials close to the Mott insulator/Metal phase boundary. 

C. Origins of Exhibited Phases 

Along the temperature-descending direction, we summarize our understanding by the 
following statements: 

(1) The staggered Q x order is a consequence of the cooperative JT effect, i.e. a combined 
effect from local JT interaction and lattice elastic energy. 

(2) The energy difference between (jr, n, it) Q x and (tt, 7r, 0) Q x orders is very small, of 
the order of meV. 

(3) The uniform Q z order is a consequence of the staggered Q x order, arising from the 
cubic term in lattice energy. 

(4) The uniform — (+)Q Z order reduces (enhances) the inter-layer AF coupling and de- 
creases (increases) the Neel temperature. 

To justify the first statement, we perform the calculation without cooperative Jahn- 
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Teller effect and obtain an orbital ordering temperature of 0.06t (~ 330K) which is far too 
low compared to the observation. The cooperative Jahn- Teller coupling arising from the 
corner-shared octahedra facilitates the staggered Q x order. The physical picture is quite 
straightforward - a. Q x distortion on one site induces an — Q x distortion on neighboring sites 
in the same x — y plane. 

The 2nd statement concerns the energy difference between the (tt, n, 0) and the (jr, n, tt) 
Q x order. There are two possible sources. The first one is the structure in the lattice 
contribution. The simple form of the lattice Hamiltonian [? ] we used in our numerical 
calculation has the same restoring force for both distortions, but a more general form given 
in Eq(j2]) will distinguish them. Assessing this possibility requires a DFT calculation of 
phonon spectrum as is discussed in Ref[? ]. The second possibility is the electronic energy 
which we now estimate from the super-exchange (SE, essentially 2nd order perturbation) 
argument. The nearest neighbor (t) and second neighbor (t') hopping processes which give 
rise to superexchange are illustrated in FigQ. 

In terms of the local JT splitting of Aq, we find that the second order superexchange 
calculation yields that the energy gain for both orders are — | -j^. Therefore within SE 
approximation, nearest neighbor hopping does not lift the degeneracy. However the second 
neighbor hopping does lift the degeneracy. We find that the (tt, tt, 0) Q x state gains — 4-J^ 
more energy than the (7r,7r,7r) state. From our DFT study [? ], t' ~ 0.035eV and Aq ~ 
1.4eV, therefore the energy difference between these two orders is of the order meV which 
is very small compared with other energy scales in the problem. We therefore believe that 
the lattice effect is dominant. We model this by allowing only the (tt, tt, 0) order in our 
calculation. 

Within our approximation, the uniform Q z order is induced by the local Q x distortion 
via the cubic term[? ] in lattice energy, so the strength of the Q z order is proportional to 
Q x . Fig([7j) shows the magnitudes of the staggered Q x and uniform Q z orders at T = O.lt. 
For these parameters we found Q z ~ 0.2Q^. 

The 4th statement concerns the relation between the magnetic order and orbital order. In 
particular it is the answer to the question that in the presence of a large staggered Q x order, 
how a small uniform Q z order affects the magnetic order. We found that a small +/-Q Z 
order can change the Neel temperature Taaf by as much as a factor of two. This effect can 
be qualitatively understood by comparing the effective magnetic couplings J, (i = z out- 
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FIG. 6: Illustrations of superexchange processes in presence of different (jr, tt, 0) and (7r,7r,7r) Q x 
orbital order. A and B are two sublattices occupying local orbitals (|1) + \2))/y/2 and (|1) - |2))/\/2 
respectively. The hopping matrices are directional and the explicit forms are given in Ref[? ]. 
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FIG. 7: The staggered Q x (solid curve) and the uniform Q z (dashed curve) orders as a function 
of doping at T = O.lt. In the bulk, the uniform —Q z order is induced by the cubic term in lattice 
energy and the magnitude is proportional to Q\. For the chosen parameters here Q z ~ 0.2Q^ 
(dotted curve). 

of-plane and x in-plane) for different orbital orders using super-exchange arguments. The 
starting point is that for each site the electron occupies the orbital \9) = cos#|l) + sinf9|2) 
(0 < 6 < 7r) which is the ground state of — (Q z f z + Q x f x ). Since the strength of Q x order 
is at least 3 times larger than that of Q z order (see Section III G), we consider Q Z /\Q X \ 
ranging from -0.3 to 0.3. As shown in FigflS]) in the presence of the staggered Q x order, the 
system is divided into two sublattices A and B on which the electron occupies orbital \9 a) 
and \9 B ). Denning cos 29 = Q z /^/Q 2 x + Q 2 z , sin 29 = Q x /y/Q 2 x + Q 2 z , one finds the occupied 
orbitals at A and B are \9), | —9). Using the 2nd order perturbation, one estimates the 
magnetic couplings from the energy difference between FM and AF spin configurations 
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FIG. 8: Illustration of the (ir, w, 0) Q x and a uniform Q 2 orbital order. The local orbital states are 
same in z direction (\6a) or \9b)) while are alternate between \9a) and \9b) in the x — y plane. 

(J = £?fm — Eaf) as 

cos 2 6 1 sin 2 6 cos 4 6 1 cos 2 # sin 2 6 1 
= i + ^ + 



T 3/16 + cos 2 sin 2 6 cos 2 0-3sin 2 # 3/16 + cos 2 9 sin 2 , , 

U9) = A + ~JK + A TA (18) 

where A jt and Auund are orbital and magnetic splitting respectively. From previous discus- 
sion we found Ajt ~ 2eV~ At, Ajiund ~ 4eV~ 81 The corresponding results are given in 
Fig© where 9/n = 0.2, 0.25, 0.3 correspond to Q Z /\Q X \ = +0.3, 0, -0.3 respectively. We see 
that for these 9 values the in-plane magnetic coupling is always FM while the out-of-plane 
changes from FM to AF when 9/ir varies from 0.3 to 0.2 (zero coupling at 9 = 0.227r). This 
SE estimate therefore implies a positive Q z order is required to produce the observed AAF 
order. In our DMFT calculation, we always find the AAF order at low temperature, but we 
indeed find the Neel temperature drastically (50%) increases when we go from small — Q z 
to small +Q Z order. Thus the trend of variation of Tjy with strain is correctly captured 
by the superexchange calculation, but other processes also contribute the overall sign. This 
result indicates that the interaction is not strong enough to justify the superexchange ap- 
proximation but that the superexchange results does capture one aspect of the important 
physics. 
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— Jx — in-plane coupling 
Jz — out-of-plane coupling 




FIG. 9: (Color online) The out-of-plane (J 2 ) and in-plane (J x ) magnetic coupling estimated from 
Eqn (ll8p with Ajt = 4t, A-Hund = 8t. The negative sign favors FM coupling. 9 = 0.25-7T represents 
the case without uniform Q z order. Positive and negative Q z orders correspond to the region 
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FIG. 10: (Color online) (a) In-plane DC resistivity p X x(T) as a function of temperature for x = 0.3. 
(b) In-plane optical conductivities for x = 0.3, T = 0.12t to 0.04t. 



VI. CMR REGIME AND CROSSOVER 



A. CMR Regime (0.3 < x < 0.5) 



We choose x = 0.3 (N = 0.7) as being representative of the CMR region. For this doping, 
the system goes from PO/PM to PO/FM around 275K. The transition is accompanied by 
an insulator/metal transition which is shown in Fig(l) of Chapter 1 in Ref[? ]. At this 
transition the dp/dT changes sign. Fig ffTOl (a) shows the calculated p{T). We indeed find 
a M/I transition across the Curie transition. This M/I transition is also reflected in the 
optical conductivities shown in Fig f[TU]) (b). Our calculations qualitatively agree with the 
experimental data of Lao.7Cao.3Mn03. To be more quantitative, Table I compares the kinetic 
energy defined as K = (ff)f Jq 7 ^ cr(u)duj obtained from both experiments [? ] and our 
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theoretical calculation. The reasonable agreement suggests that our model well captures the 
main physics (right degree of freedom and reasonable effective interactions) below 2.7eV. 
However there are several differences between calculation and data. First the experimental 
values are systematically larger. Two reasons are that the experiments inevitably involve 
transition from lower oxygen bands to Fermi level which is not included in our model and 
our calculation yields orbital order at low T which does not occur in the actual material. In 
addition our calculation overestimates the density of states around zero frequency at high 
temperature PM phase. In terms of DC resistivity, it means the high-T PM phase is not 
insulating enough. We shall discuss the possible physics accounting for this inconsistency 
later. 





Expt (Ref[? ]) 


Calculation 


FM 


0.22 


0.152 (T = 0.04*) 


PM 


0.1 


0.076 (T = O.lt) 



Table I: Kinetic energy in the unit of eV obtained from both experiments and our 

calculation, using t = 0.5eV. 

One issue from earlier calculations is that the T c for pure DE model is roughly 3 times 
higher than the observed one[? ]. The Curie temperature T c with Jaf — obtained here 
(roughly 0.08*) is ~ 40% lower than that of [? ]. Introducing G-type AF coupling Jaf(~ 
0.01* from our fit) further reduces T c to 0.075* ~ 412K, not too far from the experimental 
value ~ 275K. We expect that a large fraction of the remaining difference arises from spatial 
and thermal fluctuation effects not captured by our mean field theory. 



B. The Crossover Regime (0.1 < x < 0.3) 

In this doping range, when the temperature is lowered, the system goes from PO/PM 
to OO/PM, then to OO/FM phase and we take x = 0.2 (N = 0.8) as a representative 
doping. In this region, both cooperative JT and DE mechanisms are important. These two 
mechanisms are competing and not compatible in the following sense - the cooperative JT 
tends to break the in-plane symmetry which facilitates the staggered Q x order and localizes 
electrons, while the DE wants the system to be uniform and delocalizes electrons. This 
competition is shown in Fig([TT])(a) where the magnitude of the staggered Q x order is plotted. 
We find that when lowering the temperature, the magnitude of Q x order increases above 
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the Curie temperature, and then quickly saturate below T c . If we force the PM solution at 
low temperature, then the staggered Q x order keeps on increasing as T is decreased. This is 
consistent with the pair distribution function (PDF) measurements [? ] which show below 
the Curie temperature at x — 0.25, the peak associated with JT distortion decreases when 
lowering the temperature. Fig (ITT|) (b) shows the resistivity as a function of temperature. We 
see that the system is an insulator at high temperature and a downturn in p(T) happens 
at the Curie temperature, below which DE effects gradually takes over and the system is 
metallic. Finally we point out that around x = 0.3, T QO and T c happen around the same 
temperature (around O.lt in Fig(jIJ). We do not resolve the behavior around this point 
carefully 




0.14 



Temperature (t) Temperature (t) 

(a) (b) 

FIG. 11: (a) The amplitude of staggered Q x order as a function of temperature. The dashed 
curve are calculated at paramagnetic phase, (b) Resistivity for x = 0.2 (N = 0.8) as a function of 
temperature. The vertical dashed lines indicate the transition temperatures. 

VII. DISCUSSION 
A. Summary 

With a fixed set of parameters, our calculations semiquantitatively produce the ob- 
served phase diagram - the relative positions of magnetic and orbital orders in the doping- 
temperature plane are consistent with experiments. In particular the magnetic transition 
temperatures (both Neel and Curie temperatures) are in reasonable agreement with data 
with calculated values being about 1.5 times higher than the measured values. Some part 
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of the difference arises from the fluctuation corrections to the mean field theory, which are 
typically of the order of 30% in three dimension. As for the excitation spectra, our results 
are consistent with observed optical conductivity. In particular we reproduce the peak po- 
sitions (this is how we fit some of the parameters) and the corresponding amplitudes for a 
wide range of doping and temperature. We believe these agreements to experiments indicate 
that our model and fittings capture the essential physics of the manganite problem. In this 
section we give a more detailed discussion on several issues and on inconsistencies to data 
regarding to our results. 



B. Role and Effect of GdFe0 3 Rotation 



Our results indicate the local interaction strength is only slightly stronger than the critical 
value of Mott transition [? | implying the system is very sensitive the hopping t. As 
discussed in Ref[? ? ], the hopping is very sensitive to the structure. In particular the 
manganites form in a distorted version of the ideal perovskite structure. The most important 
important distortion appears to be a GdFeC>3-type rotation which buckles the Mn-O-Mn 
bond. Table II summarizes the relation between the bond angle, the cation composition, 
and the hopping. For a perovskite material AMnC>3, the Mn-O-Mn bond-angle as a function 
of A-site composition is taken from Ref [? ] and the corresponding hoppings are calculated 
in Ref[? ]. From this table we infer that using the same hopping t = 0.5eV for LaMn03 and 
Ca doped manganite is reasonable, but is not for Sr-doped. 



A-site 


bond-angle 


hopping (ratio) 


Lax Ideal 


180 


0.65eV (1) 


Lao.7Sro.3 


166 


0.58eV (0.9) 


Lao.7Cao.3 


160 


0.53eV (0.81) 


Lai Real 


155 


0.5eV (0.78) 



Table II: The composition of A-site elements, its corresponding Mn-0-Mn angle and 

effective hopping t. 

In Ref[? ], we show that for t = 0.65eV LaMn03 is not a Mott insulator. When using 
t = 0.6eV to simulate Lao.7Sro.3Mn03, we find that (1) the Curie temperature increases 
from ~ 420K (t = 0.5eV) to ~ 530K (shown in the dashed curve in Fig(jl])) and (2) the high 
temperature PO/PM phase becomes almost metallic (dp/dT is very flat, and the minimum 
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FIG. 12: The optical conductivities calculated for parameters which simulate Lao.7Sro.3Mn03. 
Solid and dashed curves are computed above and below the Curie temperature. To convert the 
frequency into physics units [eV], divide the x-axis by two. 

around a(u = 0) almost vanishes). The optical conductivities for t = 0.6eV for two temper- 
atures, just above and below T c , are shown in Fig f[P2l) . Both of our findings (T c and <j(uj)) 
are consistent with the difference between Lao.7Sr .3Mn03 and Lao.7Ca .3Mn03 reported in 
Ref[? ]. We emphasize, however, that the main message here is that for the given local 
interaction strength, the system is very sensitive to the bandwidth and any uncertainty in 
estimating parameters could easily drive the system to either Mott insulating or metallic 
phases. 

C. Orbital Ordering 

We found that with our standard parameters, the calculated T OQ decreases too slowly 
as a function of doping x (see Fig(Jl|)). Fine tuning parameters (e.g. varying e and Uq) 
can correct this problem but this degree of data fitting is somewhat arbitrary so we do not 
persue it here. 

However, we point out that using our standard parameters, the orbitally disordered and 
orbital ordered phases behave very similarly as far as the excitation is concerned. To demon- 
strate this we show in Fig(fT3l) the optical conductivities for x = 0.3, T = 0.08t at both OO 
and PO phases. One sees that the difference is very limited. Furthermore we found the 
Curie temperature obtained from these two phases are very close (difference < 5%). These 
results suggest that the crucial electronic physics is controlled by local distortions which 
(because the correlations are strong) are well formed. The important effect caused by the 
orbital order is actually the uniform Q z order which introduces an isotropy to the system 
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FIG. 13: (Color online) The optical conductivities calculated at x = 0.3, T = 0.08t for PO (solid) 
and 00 (dashed) phases. 



and whose sign substantially affects the Neel temperature for the undoped case. 

We also mention here that the effective Mn-Mn hopping being through oxygen p orbitals 
also introduces an intersite orbital coupling [? ] which is very similar to the cooperative 
JT effect and is referred to as the "charge-transfer" mechanism. This can be understood 
from the super-exchange argument where we consider a simple Mn-O-Mn system and com- 
pare energies of different orbital configurations by the perturbation expansion of the Mn-0 
hopping t P d [? ]. In the model where the oxygen orbitals are not included, those "virtual 
processes" lead to an spin-independent orbital-exchange interaction [? ] as 

H „ ex = A I? 1 ?-,* (19) 

i,ct=x,y,z 

with I z = t z , I x = —\t z — ^t x , P = —\t z + ^t x , and A a positive coefficient. The simple 
mean field approximation to the term produces an external field on site o as 



H _ ex = a io((n+a+(n- 



a=x,y,z 

A 



1 l 

T Z -^;(2E Z — E x — E y ) + t x -j=(—E x + E y 



(20) 



with E z = ^J2 a =±z( T z) a and E x{y) = ^E Q=±I (± v )[^)a + (-)V3(r x ) a }. which is of the 
same form of the cooperative JT effect derived in the appendix. Therefore in our approxi- 
mation where the orbital order and structural JT distortions are equivalent, including the 
charge-transfer mechanism amounts to a reinterpretation of our cooperative JT parameter 
e = 2A and does not change any results. 
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D. High-T Insulating Phase 

Using the standard parameters, our calculation obtains an insulating behavior at high 
T PO/PM phase for doping ranging from x = to x ~ 0.4. With the semiclassical ap- 
proximation, the electron-electron interaction is replaced by some classical fields and the 
impurity problem becomes polaron-like [? ? ]. The high T insulating phase away from 
zero doping should be therefore interpreted as a phase separation between N = 1 orbitally 
fully-polarized state and N = state. Since our estimate indicates [? ] the on-site Coulomb 
interaction is roughly three times stronger than the electron-lattice, the semiclassical method 
may overestimate the insulating behavior under single-site DMFT approximation. 

The other issue is that compared to the experiments, our calculations overestimate the 
optical conductivity around zero frequency at high temperature PO/PM phase. This might 
be due to short-ranged correlations not included in the single site DMFT approximation. 
According to recent cluster DMFT studies of the 1-band Hubbard model [? ? ], including 
the short-range correlation significantly reduces the low energy density of states. We also 
observe that in doped CMR systems there is a strong empirical association between insulat- 
ing behavior (dp/dT < with p DC resistivity and T temperature) and strong short-ranged 
Jahn- Teller (polaron glass) order [? ? ? ]. Including spatial correlations beyond the single 
site approximation is an important topic for future regard. 

E. Missing Phases 

Our calculation misses two phases. First around x = 0.5, a charge ordered (CO) phase 
occurs, accompanied by one particular orbital and magnetic order called CE phase [? j 
which requires a very large unit cell (4 x 4) in the x — y plane. Our in-plane unit cell is 
not large enough to include this phase. However at x = 0.5 we do find that the convergence 
becomes more and more difficult when lowering the temperature (below T = 0.04t) which 
may be an indication of CE phase. Second, around x — 0.1 — 0.2 we do not get the FM 
insulating phase at low temperature. 
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F. Limits of Approximation 



Now we discuss the limits of our approximation. First we discuss the breathing mode 
polaron effect. In the current approximation the breathing mode coupling is treated in 
simple mean field and therefore has no effect in the charge- uniform phase. To include the 
breathing mode polaron, one has to consider the fluctuation of the breathing-mode distortion 
by integrating over Q Q field when computing the impurity model [? ]. Since the real 
time-consuming computation involved in our approach is doing multi-dimensional integral 
(see Section III.F), performing an additional integral is now beyond our computational 
power. It is possible that the breathing mode polaron is also crucial for the charge order 
at x — 0.5. Since with breathing mode polarons electrons are already localized but just 
randomly distributed at high temperature (therefore the system is not charge-uniform, but 
can be treated within single-site DMFT [? ]), the CO is then formed at low temperature to 
gain more energy from the gap. Without the polaron to localize electrons, it is very hard to 
get CO (usually it requires some ti ™ +iii ' r ; " u^^a rfr„^„™ which is not the case here). 
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FIG. 14: Spectral function calculated for one-band Hubbard model for bandwidth 6t, on-site 
U = 16t at T = O.lt. For both half filling (solid) and N = 0.8 (dashed), the lower and upper bands 
have the equal weights. 

Finally we point out that the SCA does not treat the quantum physics of the Mott 
insulator faithfully. To be specific, we take 1-band Hubbard model with strong coupling 
as an example. With the SCA the metal/insulator transition occurs at U ~ U c \ so the 
effect pf the Kondo peak is absent. Further, independent of doping x, the upper and lower 
Hubbard bands have the same weight with SCA while in the reality, the upper Hubbard 
band represents adding one electron to the occupied site the whose weight is ought to be 
1-x. This consideration implies the SCA works well at half filling and becomes less reliable 
away from it. This is illustrated in Fig (fl~4|) which show the spectral functions for N = 1 
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(half-filling) and N = 0.8. We see that in both cases the upper and lower bands have 
the same weight. For our 2-band manganite model in the strong coupling limit, the SCA 
solution for PO/PM phase results in 4 bands with weights 1-x, 1+x, 1+x, 1-x (from low 
energy to high) respectively. However the 3rd peak corresponds to adding one electron to 
the state with same orbital but opposite spin whose weight should be 1-x. Based on the 
same argument we conclude that the SCA for 2-band model is more reliable without doping 
than with doping. A more accurate treatment of the doped phase requires an improved, 
fully quantum impurity solver. 

VIII. CONCLUSION 

A general model for bulk manganite, including electron-electron, electron-phonon, and 
phonon-phonon interactions is formulated and solved by semiclassical approximation. Our 
calculation is qualitatively good in the sense that it yields the right distribution of phases in 
the (x, T) plane and produces the correct low energy excitations as described in section 2.5.2. 
The physical origin of each exhibited phase is identified within our model. For the LaMn0 3 
below T OQ , the exhibited in-plane staggered Q x order is mostly driven by the cooperative 
Jahn- Teller (lattice effect) rather than the pure electronic effect while the uniform — Q z order 
is a consequence of the anharmonic term in lattice energy. Our results indicate that the local 
interaction strength is only slightly stronger than the critical value for Mott transition and 
the system is consequently very sensitive to mechanisms controlling the effective bandwidth. 
With this local interaction strength, the orbitally ordered and orbitally disordered phases 
behave very similarly. As the doping increases, the electrons start to delocalize and after 
x > 0.3, the double-exchange mechanism dominates so orbital order vanishes and the system 
has the FM/Metallic ground state. 

Discrepancies between our calculation and the observations are also carefully discussed. 
In particular our calculations overestimate the optical conductivity around zero frequency 
at high-temperature insulating phase. This inconsistency leads us to conclude that the key 
physics we are missing in the calculation is the short-ranged correlation. In the future we 
will include the short-ranged correlation and also adopt a better impurity solver for this 
problem. 
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APPENDIX A: THE VALIDITY OF SEMICLASSICAL APPROXIMATION 

In this appendix we examine the validity of the semiclassical approximation by comparing 
the excitation spectrum computed using the SCA results to the exact eigenstates of the local 
Hamiltonian. As discussed in Section II, we assume that the crystal field (ligand field) is 
large enough that the t 2g levels are in their maximum spin state and that the pair hopping 
between t 2g and e g orbitals is quenched. In this case the on-site Hamiltonian in the e g 
manifold is 

0,0' i=l,2 

- 2Jsx • s 2 - 2J H S C • (si + s 2 ) + A(rai - n 2 ) (Al) 

Here Sj = ^ a /3 c L^a/3 c ^> $c nas magnitude 3/2 and A is the crystal field splitting arising 
from the long-range Jahn- Teller order. In spherical symmetry J H = J; we assume this 
henceforth. The eigenstates are characterized by the particle number, total spin and total 
e g spin, and the orbital configuration. There are 16 1-electron and 24 2-electron eigenstates, 
taking the configurations of the core spin into account. 

To compare the exact solution of the local Hamiltonian to experiment and the semiclas- 
sical calculation, we need the quantity AE(S) = E(n = 2, S) + E(n = 0, S = 3/2)-2E(n = 
1, S — 2) which gives the locations of peaks in the optical conductivity in the atomic limit. 
The following table lists the eigenstates and the corresponding transition energies. 



States 


AE 


Semiclassical 


M 2 (5/2) (6) 
M 2 (3/2) (4) 
M 2 (l/2) (2) 


U - 3 J/2 + 2A 
U + 7J/2 + 2A 
U + 13 J/2 + 2A 


2(U S + J H ) + 2U Q 
not accessible 


^-(3/2) (4) 
M(3/2) (4) 
1 E + (3/2) (4) 


U + 9 J/2 + 2A - V4A 2 + J 2 


2{U S + J H ) 
not accessible 
2(U S + J H ) + 2U Q 


U + 9 J/2 + 2A + V4A 2 + J 2 
U + 7J/2 + 2A 
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Table III: The 2-electron eigenstates and 
the corresponding transition energies 



Determining the coupling strength in Eq ljAlft by fitting the optical data [? ] is described 
in detail in Ref[? ]. Here we simply quote the results, U = 2.3 ± 0.3eV, 2 A ~ J ~ 0.5eV. 
Following the analysis and notations in Ref[? ], there are three optical peaks located at 

AE HS = U-3J/2 + 2A 

AE LS = U + 9 J/2 + 2A - V4A 2 + J 2 

AE+ S = [7 + 7J/2 + 2A (A2) 

We now compare this result to the semiclassical calculation. From Fig(j3j) we observe 
three peaks in the optical conductivity: a low-lying peak at energy 2Uq which we identify 
with AEhs, an intermediate peak at energy U s + Jh which we identify with AEfg, and a 
higher peak at the sum of these energies. This highest peak represents physically the states 
1 E + (3/2) and 3 A 2 (3/2) where both orbitals are occupied while the total spin (including the 
core spin) is 3/2. The 3 A(3/2) state in large A limit represents a state where both electrons 
occupy energy-disfavored orbital which cannot be reached by a single hopping and has no 
correspondence in the SCA. It is the defect of the semiclassical approximation that the 
highest peak is too high in energy. However this defect is not serious because the high-lying 
states are not important for our analysis. 

APPENDIX B: EFFECTIVE POTENTIAL 

In this appendix, we describe in detail how we encode the inter-site lattice coupling into 
the single-site impurity problem. The basic logic is the following. First we write down the 
energy functional for the lattice problem in terms of fields labeled by site index 0j which 
couples to some local quantity pi, then the local partition function is obtained by integrating 
out all fields except the field at origin site 0o- The long-range order corresponds to some 
spatial pattern of pi which generates an extra coupling to local field 4>q. This extra coupling 
depends on the long-range order containing information from other sites pi i ^ o. We first 
give a general functional for lattice elastic energy then work out 1-dimensional case explicitly 
with a specific lattice model. Finally we derive the formalism used in our calculation. 
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1. General Functional of Elastic Energy 



The goal here is to derive the elastic energy in terms of three even-parity MnOg distortion 
modes. As mentioned in the text, the lattice degree of freedom includes oxygen motion 
along Mn-0 bond and manganese general displacement 8^. Assuming the spring constant 
between adjacent Mn-0 is l/Ki, a general elastic energy is 

E lat = -L^ 5 t-<) 2 + ( s t ■■" ^ 



u i-u 



2K X 



v 7 fc ~ k 2 

k,ab k,ab 



S ' / L k —k k —k y ' k —k 1 



k.a 



\ E«\k)5l8K + \ £ D°»{k)ulu»_ % (Bl) 



K ' k ~ k 2 . 

k,ab k,ab 



where E ab (k), D ab (k) represent general harmonic coupling Mn-Mn, 0-0 displacements, and 
a, h sums over x, y, z. To get rid of the Mn motions, we use the saddle point approximation 



= which leads to 

— k 



%=^Y)l + MZ)/2)]*\l + ^>i (B2) 

b 

and the lattice energy in this approximation is Ei at = Ylk a b U k^ 1 

(k)u b - with 

fi ab 1 -> 1 

rh ab (k) = — - + e-^)[l + E(k)m]-J(l + e ik >) + -D ab (k) (B3) 



Defining strain variables vf = u\ — u\_ a , vi = — e we express E iat in terms of v| 
which is 



E lat = J2v a % m ab (k)v b _ % (B4) 

k,ab 

where m ab (k) = 1 _ f ^ ika fh ab {k) ■ The advantage of expressing E\ at in strain variables is 
that they are closer to the even-parity distortion modes defined in Eqn([S]). m ab (k) allows us 
to estimate the proximity effect for structural order. In particular if we are interested in how 
(n, 7r, 0) Q x order propagates along z direction, then the quantity to study is m ab (ii, tt, k z ). 
The explicit form of m ab (k) is model-dependent and here we only consider spring constants 
between adjacent Mn-0 (l/Ki) and Mn-Mn (I/K2) which are of most importance. 
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2. 1 Dimensional Mn-O Chain 



Now we explicitly work out the local effective potential in the 1-dimensional case. The 
procedure is outlined here. We first adopt the procedure described in the previous subsection 
to express (with saddle point approximation) the elastic energy in terms of strain variables 
in real space V{. Then the effective potential is obtained by integrating out all Vi except the 
one at origin vq. For ID Mn-0 chain, we drop the index a, b since there is only one direction 
and the elastic lattice energy is 

E lat = JL - u % f + (5, - u^f] + JL J2(5 l+1 - 5 t f 

1 i 2 i 

1 . 2 . 

= -rr ^y\u k u_ k + 5 k 5- k - u k 6- k (l + e ik )] + — V sm 2 (k/2)S k 6. k (B5) 

Kl k K2 k 

For this case, -^rp = sin 2 (/c/2) and Dik) = 0. The saddle point approximation = 
implies 

u k (l + e* k ) 

4 " 2 + 4^sin 2 (A;/2) (B6) 
where K = Ki/K 2 . The effective energy functional E iat (as a function of u k only) is therefore 



2K+l^ r sin 2 (A;/2) 
L "" ~lcr^l + 2Ksm 2 (k/2) Ukl1 -'' 

2K+1^ 1 



^ l + 2KL 2 (k/2) VkV - k (B7) 



4ATi ^ l + 2ATsin 2 (fc/2) 

To see how local strains at different sites couple to one another, we express E\ at in the 
real space Uj. 

2g + 1 ^ 1 
Elat ~ ~AKY V l + 2irsin 2 (A;/2)^- fc 
1 2i? + 1 g^^-r,) 



N AK 1 l + 2Ksin 2 (fc/2) 



2/7 /o\ ^J 

- — — ± -r sir 

i,j k 

2K+1 



AK X 



Y^f(i-j)vi»i (B8) 

where /(i - j) = ^ E fc 1+2 g^ /2) = ^ £, dk 1+2 ^g fc/2) (£ fe - for lattice 

constant a = 1). For this simple model, the integral can be done analytically (the most 
straightforward way may be changing variable z = e tk and then using the residue theorem!). 
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By defining a = 1 + K 2 / K\ , E iat becomes 



hi 



= ^E e '^ ( B9 ) 

where if = 1K\j\/lK + 1, and e = a — V« 2 — 1 < 1- Note that the coupling between local 
strains is exponentially decay since e n = e~ a ' n ' with a = — lne. 

Including the the electron-lattice coupling J2 i hi v ii the total energy is 

E = E lat + E e ^ = A H v i v j + hiVi ( B1 °) 

ij i 

where hi in this case is the charge density at site i. The effective potential at site is given 
by integrating out the degrees of freedom of all other sites Vi,v 2 , ...,i>jv, i.e. 

/ / / 



J dv e- v °" M = J dv e- AooV o Jd Vl ..dv N exp[-J2 A oiVoVi-J2 Ai ^~J2 hi 

r i ' i ' i ' 

= J dv exp[-(A 00 - - A AoiA^Ao^vl + - ^ M^A)^o + h i A ^ h i 

where Yl' means site is excluded in the summation. 
The effective potential is 

I ' 

V eff (v ) = Dvl - - hjA'/AoiVo + const (Bll) 

where D = A 00 — A 0i A ij~ A 0j- We see that the charge density at site hi(i ^ 0) also 
contribute to the "external" field coupling to v . 



3. 3 Dimensional Case 

The ID result can be easily generalized to the 3D case. For the model we considered, the 
lattice energy in k space is 

2K + 1 1 



lot 



^I^H^tfiw^ (B12) 



From this expression, we find that in our simple model there is no proximity effect for (it, n, 0) 

— * 

order of any kind since there is no coupling between different components of k. We also 
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notice that the energy cost is at its minimum when k a = tt (staggered order of any kind), 
therefore at integer occupancy the system prefer some staggered long-range order since the 
staggered order lowers electronic energy. 

Assuming e is small thus only including the nearest neighbor coupling, the lattice energy 
in real space is 



Eiat = ^ R £ Mf + 2ev?v? +a ] 



(B13) 



i,a=x,y,z 

One can also express E[ at in three Mn0 6 even parity modes Q by the following transformation 



to ^ 

Vi,0 
Qi,x 

y Qi,z J 



( i 
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JL L n 

1_ 1_ _2_ 
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i,x 



v i,y 

\ V i,z ) 



However it is more convenient to work in stain field v until we obtain the local effective 
potential which will be expressed in Q. 

The electronic source fields hi are defined as 

hi, = (eab4,a C i,b) 
K,x = ( T ab C i,a C i,b) 
Kz = (Tfcl Ci th ) 



and the local electron-lattice coupling is 



Ej T = -h Qo - (h x Q x + h z Q z ) = -(v x H x + v y H y + v z H z ) 



(B14) 



with 



H x = ±h + -±=h x - ±h z 
h z = + j=h, 



(B15) 



Following the procedure for 1-dimensional case, the effective potential at the origin site 
is therefore 

1 



Veff = D(v 2 x + v 2 y + vi) + -E-v 
= D(Ql + Ql + Ql)+ € -F.Q 



(B16) 
(B17) 
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with E xM = Yl H U y , z )AyMiX,z)My,z) A °^y>z) and ^ = u ^- Vi/,*) labels the sites alon S 
x{y, z) axis. We call F effective external field. Keeping only the linear term in e, we have 

B= _L (1 _ £7 4)~^ (B18) 

E x H+ xx -\- H_ xx 

E z = H +Z ^ z + H_ zz (B19) 

(B20) 

and F = (j=(E x + E y + £ 2 ), j=(E x - E y ), ^(-E x - E y + 2£*)). 
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